Transition from the Couette- Taylor system to the plane Couette system 
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We discuss the flow between concentric rotating cylinders in the limit of large radii where the 
system approaches plane Couette flow. We discuss how in this limit the linear instability that leads 
to the formation of Taylor vortices is lost and how the character of the transition approaches that of 
planar shear flows. In particular, a parameter regime is identified where fractal distributions of life 
times and spatiotemporal intermittency occur. Experiments in this regime should allow to study 
the characteristics of shear flow turbulence in a closed flow geometry. 
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The transition to turbulence for a fluid between con- 
centric rotating cylinders has attracted much experimen- 
tal and theoretical attention. Ever since Taylor's success 
|E| in predicting and observing the instabilities for the 
formation of vortices the system has become one of the 
paradigmatic examples for the transition to turbulence 
and a large number of bifurcations have been analyzed 
in considerable detail ^ ||. The limiting case of large 
radii and fixed gap width where the effects due to cur- 
vature become less important and where the system ap- 
proaches plane Couette flow between parallel walls has 
received much less attention. In this limit the character 
of the flow changes: plane Couette flow is linearly stable 
and the mechanisms that drive the transition to turbu- 
lence are still unclear. The question we address here is 
to what extend the Couette- Taylor system can be used 
to gain insight into the dynamics of plane Couette flow. 

This problem is of both experimental and theoretical 
interest. As mentioned, the experimental situation for 
Couette- Taylor flow is much better, there being numer- 
ous facilities and detailed studies of patterns, boundary 
effects and critical parameters . The moving bound- 
aries in plane Couette flow reduce the experimental ac- 
cessibility and the possibilities of applying controlled per- 
turbations. On the theoretical side it is an intriguing 
question how the change in stability behaviour from the 
Couette Taylor system to the plane Couette system oc- 
curs. Studies by Nagata |Q show that some states from 
the rotating plane Couette system survive the limiting 
process and appear in finite amplitude saddle node bifur- 
cations in the plane Couette system (see also the investi- 
gation of this state by Busse and Clever S). Unless the 
transition from linear instability dominated behaviour in 
Couette- Taylor flow to the shear flow type transition in 
plane Couette flow is singularly connected to the absence 
of any curvature it can be expected to happen at a finite 
radius ratio near which interesting dynamical behaviour 
should occur. 

We should mention that there are other useful embed- 
dings of plane Couette. Busse and Clever || start from a 
layer of fluid heated from below with cross flow and pro- 
ceed to study the stability and parameter dependence of 
the states. And Cherhabili and Ehrenstein |j) start from 
plane Poisseuille flow and find localized solutions, albeit 



at Reynolds numbers higher than the ones studied here. 

Our aim here is to follow some of the instabilities in the 
Couette- Taylor system to the limit of the plane Couette 
system and to identify the parameters where the change 
in behaviour occurs. In particular, we study the transi- 
tion from laminar Couette flow to Taylor vortices and the 
instability of vortices to the formation of wavy vortices. 
Note that the asymptotic situation of plane Couette flow 
can be characterized by a single parameter, a Reynolds 
number based on the velocity difference, whereas Couette 
Taylor flow has at least two parameters, the Reynolds 
numbers based on the velocities of the cylinders. This 
extra degree of freedom provides an additional parame- 
ter that can be used to modify the flow without changing 
the basic features. 

In cylindrical coordinates (r, (f>, z) the equations of mo- 
tion for the velocity components (u r , u$,u z ) can be writ- 
ten as 
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where the modified Nabla and Laplace operators are 
V = e r d r + ej-d,h + e z <9 



A = d rr + ^-<9<>0 + d zz , 



(1) 

(2) 

(3) 
(4) 

(5) 
(6) 



and where e, are the unit basis vectors |10| 

The terms in eqs. ((j])-(f|) are arranged so that all the 
ones on the right hand side vanish when the system ap- 
proaches the plane Couette system, i.e. in the limit of 
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large radii but finite velocities at the cylinders. The re- 
maining ones become the equations of motion for plane 
Couette flow in cartesian coordinates (x, y, z) if the iden- 
tification x — r and y = (fyr is made. However, there are 
other ways of taking the limit of a small gap that lead 
to different limiting systems. For instance, the case of 
almost corotating cylinders with high mean rotation rate 
gives rise to plane Couette flow with an additional Corio- 
lis term ('rotating plane Couette flow' j7J). Another limit 
corresponds to the case of counterrotating cylinders with 
diverging rotation rates pd[ |. In our numerical work we 
use the full equations, without any reduction in terms. 
This allows us to extend Nagata's work from the rotating 
plane Couette flow to the full Couette- Taylor system. 

The velocities at the inner and outer cylinder (distin- 
guished by indices i and o, respectively) are prescribed 
and define the boundary conditions 



u<j,{r = R x ) = tt x R x , 

u r (r = R x ) = u z (r = R x ) = 0, 



i. o . 



(7) 
(8) 



For the choice of dimensionless quantities we appeal 
to the plane Couette flow limit. There the relevant 
quantities are the velocity difference between the walls, 
AU = Ri^li — R ^i , and the gap width d = R — R4. 
Without loss of generality we can always assume > 0. 
The Reynolds number for plane Couette flow is based on 
half the velocity difference and half the gap width, 



Re = 



AUd 



(9) 



For the Couette Taylor system there are two Reynolds 
numbers based on the gap width and the rotation rates 
of the inner and outer cylinders, 



Re x — R x £l x d/v , 



(10) 



where the index x can stands for i or o, the inner and 
outer cylinders. The plane Couette flow Reynolds num- 
ber thus is Re = (iie, — Re )/A. The ratio of these 
Reynolds numbers will be called 



jj, = Re /Rei 



(11) 



(the tilde is used to distinguish it from /i = o /0j, a 
frequently defined quantity not used here). 



rj = RijR 



(12) 



denotes that ratio of radii. 

Experiments and numerical simulations show that 
plane Couette flow undergoes a subcritical transition to 
turbulence around Re PC F ~ 320 The Couette- 

Taylor system shows a first linear instability to the for- 
mation of vortices (Taylor- vortex flow, TVF) at Reynolds 
numbers that depend on the rotation rates and the curva- 
ture of the cylinders. In order to see shear flow dominated 
dynamics the critical Reynolds number for the linear in- 
stability has to be above RepcF- The formation of TVF 



occurs at Reynolds numbers that can be parametrized in 
the form 



Re = A(H)(l-r}) 



-1/2 



B(/2) 



(13) 



for 77 < 1 Jl5| |. This number is larger than the tran- 
sitional Reynolds number for plane Couette flow if 77 is 
sufficiently close to one. The minimal radius ratio 77320 
where the linear instability occurs for Re > 320 strongly 
depends on the ratio of the Reynolds numbers of inner 
and outer cylinder. A few examples for minimal radius 
ratios 77320 are summarized in Table |. 

Very important for the transition to turbulence in lin- 
early stable systems are nonlinear processes that could 
give rise to some finite amplitude states, perhaps station- 
ary or periodic, around which the turbulent state could 
form. One candidate that could serve as a nucleus for tur- 
bulence in plane Couette flow is the stationary state first 
calculated by Nagata . He observed that the wavy vor- 
tices that form in a secondary instability from the TVF in 
the rotating plane Couette system can be followed to the 
limit of the plane Couette flow where they become part 
of a saddle node bifurcation at finite Reynolds numbers. 
This state was also identified and studied in a different 
limiting process by Busse and Clever ||. They found 
that the critical axial and azimuthal wavelengths for this 
state are 



A 2 = 7r and Xa, = 2n . 



(14) 



This is roughly twice the critical wave lengths that would 
be expected for Taylor vortices. 

We developed a numerical code for the solution of the 
full Navier-Stokes equation using Fourier modes in ax- 
ial and azimuthal direction and Legendre collocation in 
the radial direction. The pressure terms were treated 
by a Lagrange method. The period in z and 4> was de- 
termined by the fundamental wave lengths ( |T4| ) of wavy 
vortex flow. 

The continuation of the wavy vortex flow from the 
Couette- Taylor system to the plane Couette system is 
shown in Fig. [I] for the case of the outer cylinder at rest 
(jl = 0) and for counterrotating cylinders with p, = —1. 
For small 77 the wavy vortex develops from a secondary 
bifurcation of TVF, but for sufficiently large 77 the wavy 
vortex state is created first in a saddle node bifurcation. 
The critical Reynolds number for the formation of Tay- 
lor vortices diverges as 77 approaches one, but the one 
for the formation of wavy vortices approaches a finite 
value. Thus the gap in Reynolds numbers between the 
two transitions widens and the region where plane Cou- 
ette flow like behaviour can be expected increases with 
77 approaching one. The radius ratios ?7 C and Reynolds 
numbers Re c of the codimension two point where the in- 
stabilities for TVF and wavy vortex flow cross are listed 
in Table |[ The ratio of radii 77 c where the linear insta- 
bility of Couette flow and of the Taylor vortex flow cross 
is a non-monotonic function of the ratio jl of rotation 
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speed. Both the critical Reynolds numbers for the linear 
instability of the Couette profile and for the formation of 
wavy vortices increase with decreasing jl, but at different 
rates and with different dependencies of rj. As a conse- 
quence there seems to be a local minimum near about 
0.93 for jj, close to —1. 

For the parameter value considered here the curvature 
of the cylinder walls is geometrically small (see Fig. ||). 
On the length of one unit cell in (^-direction the rela- 
tive displacement in radial direction from a planar wall 
is about 7r(l — rj), i.e. only 3% for 77 = 0.99. 

The critical Reynolds number for the formation of 
wavy vortex flow (WVF) seems to converge to the same 
value for both ratios jl shown in Fig. [l]. The critical 
Reynolds number as well as the rotation speed of the 
wavy vortices for several different ratios jl are collected 
in Fig. 0. The rotation speed is defined as the angu- 
lar phase velocity oj of WVF times the mean radius 
R = (Ri + i? )/2 minus the mean azimuthal velocity 
v = {uOiRi + lo R )/2. For all ratios between the speed 
of inner and outer cylinder the critical Reynolds number 
for the formation of the wavy vortex state converges to a 
value of about 125 and the speed of rotation goes to zero. 
The limiting state that is approached is the stationary 
Nagata-Busse-Clever state. The velocity field of a wavy 
vortex solution at rj — 0.993, jl = — 1 and Re = 124 is 
shown in Fig. ||; it differs little from the corresponding 
plane Couette state obtained by Busse and Clever |^J, 
both in appearance and in critical Reynolds number. 

In the region above the wavy vortex instability but 
below the linear instability the dynamics of perturba- 
tions shows the fractal life time pictures familiar from 
plane Couette flow Jig ]. Fig. ^ shows an example at a 
radius ratio of 77 = 0.993 and a Reynolds number ratio 
of jl = —1. The initial state was prepared by rescal- 
ing a WVF field obtained at very low radius ratio and 
Reynolds number. It is interesting to note that even 
with this initial condition, which is at least topologi- 
cally close to the Nagata-Busse-Clever state, it is not 
possible to realize a turbulent signal in its neighborhood: 
the state quickly leaves this region in phase space. One 
might have hoped that in spite of the linear instability 
of the Nagata-Busse-Clever state other states created out 
of secondary bifurcations could have supported some tur- 
bulent dynamics in its neighborhood, but the numerical 
experiments do not support this. The gap between the 
Reynolds number where the WVF state is formed and the 
one where typical initial conditions become turbulent is 
about the same as in plane Couette flow: the WVF states 
forms around Re = 125 and the transition to turbulence, 
based on the requirement that half of all perturbations 
induce a long living turbulent state, occurs near a value 
of Regans = 310, very much as in plane Couette flow 
& 

In summary, we have identified parameter ranges in 
the Couette- Taylor system where some of the character- 
istics of the plane Couette system can be found. These 
parameter ranges include radius ratios that can be real- 



ized experimentally. Investigations in this regime should 
be rewarding as they open up the possibility to study the 
properties of the transition in a closed geometry and to 
switch continuously between supercritical and subcritical 
transition to turbulence. The observation of a codimcn- 
sion two point where the linear instability to TVF and 
the secondary instability to wavy vortex flow cross should 
provide a starting point for further modelling of the tran- 
sition in terms of amplitude equations. 
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FIG. 1. Bifurcations to Taylor vortex flow (TVF) and wavy 
vortex flow (WVF) in Couette- Taylor flow for the outer cylin- 
der at rest (p, = 0) and counter-rotating cylinders (fi = — 1). 
The vertical line indicates the parameter range of the lifetime 
measurements of Fig. M at fi = — 1. 



FIG. 3. The convergence to the Nagata-Busse-Clever state 
for different rotation ratios fi. In the limit of r\ going to one the 
wavy vortex states for all p, approach the same flow that moves 
with the mean velocity azimuthally. The top diagram shows 
the rotation speed and the botton one the critical Reynolds 
numbers. 




2nd 




FIG. 2. Geometrical curvature of the cylinders in the Cou- 
ette- Taylor flow and the plane Couette flow limit. Shown 
is one fundamental azimuthal wavelength for different radius 
ratios r\ as indicated. 



FIG. 4. The wavy vortex flow state near the Na- 
gata-Busse-Clever state at r\ = 0.993, ft, = — 1 and Re — 124. 
Shown is only the disturbance, without the Couette profile. 
The frames from left to right show cuts through the (r, z) 
plane at azimuthal wave lengths <f> = 0, 7r/4, n/2, 37r/4 and 
7r. The vectors indicate the r and z components of the ve- 
locity field and shading the ^-component. The inner (outer) 
cylinder is located at the left (right) side of each frame. 
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FIG. 5. Life time distribution in Couette- Taylor flow at 
r\ — 0.993, p, = — 1 and for the indicated range of Reynolds 
numbers. 
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TABLE I. Parameters connected with the Couette- Taylor 
system in the limit of large radii. A and B are the coefficients 
in the parametrization ( p.3[ ) of the primary instability. 77320 
is the radius ratio where the primary instability lies above 
Re = 320; finally, r/ c and Re c are the parameter values for 
the crossing of the stability curves for Taylor vortex flow and 
wavy vortex flow. 
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